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The relativistic light-front dynamics (LFD) method has been shown to give a correct 

description of the most recent data for the deuteron monopole and quadrupole charge 

form factors obtained at the Jefferson Laboratory for elastic electron-deuteron scattering 

for six values of the squared momentum transfer between 0.66 and 1.7 (GeV/c) 2 . The 

good agreement with the data is in contrast with the results of the existing non-relativistic 

approaches. 

In this work we firstly make a complementary test of the LFD applying it to calcu- 
late another important characteristic, the nucleon momentum distribution n{q) of the 
deuteron using six invariant functions fi (i — 1, 6) instead of two (S- and D- waves) 
in the nonrelativistic case. The comparison with the y-scaling data shows the deci- 
sive role of the function /s which at q > 500 MeV/c exceeds all other /-functions (as 
well as the S- and D-waves) for the correct description of n(q) of the deuteron in the 
high-momentum region. Comparison with other calculations using S- and D-waves corre- 
sponding to various nucleon-nucleon potentials is made. Secondly, using clear indications 
that the high-momentum components of n(q) in heavier nuclei are related to those in the 
deuteron, we develop an approach within the natural orbital representation to calculate 
n(q) in (^4, Z)-nuclei on the basis of the deuteron momentum distribution. As examples, 
n(q) in 4 He, 12 C and 56 Fe are calculated and good agreement with the y-scaling data is 
obtained. 



I. INTRODUCTION 



The most recent experimental data on the deuteron structure functions and tensor polarizations ob- 
tained at Thomas Jefferson Laboratory have been reported in (and references therein). The new 
data on tensor polarization observables measured in elastic electron-deuteron scattering for the squared 
momentum transfer between 0.66 and 1.7 (GeV/c) 2 |l| make it possible to determine quite precisely the 
deuteron charge form factors for momenta comparable with the deuteron mass. Thus these data are 
important to probe the relativistic dynamics inside the deuteron. It was shown in (T| that the empirical 
results for the tensor polarization i 2 o an d both monopole (Gc) and quadrupole (Gq) charge form factors 
are in very good agreement with the calculations within two relativistic and covariant models: i) a model 
developed in the framework of the explicitly covariant version of light-front dynamics (LFD) (e.g. J7|-|l0||) 
and, ii) a model using a three-dimensional reduction of the Bethe-Salpeter equation [fj"l"| . At the same 
time the non-relativistic impulse approximation approach fl2|| (even with inclusion of meson exchange 
currents and relativistic corrections) cannot explain, for instance, the behaviour of the monopole charge 
form factor in the region of the first node and the secondary maximum shown for the first time in . 

It is well known jl3|-[l5| that one of the main feature of the realistic nuclear models (beyond the 
mean-field approximation) is that they have to describe simultaneously both important ground state 
characteristics, the nucleon momentum and density distributions of the system. This is not the case 
in the Hartree-Fock approximation where the dynamical short-range and tensor nucleon-nucleon (NN) 
correlations are only partly incorporated. The main characteristic feature of the data on the nucleon 
momentum distribution n(q) extracted by the nuclear y-scaling analysis ]l6|-|T9| (we shall call them all 
along this paper y-scaling data (YSD)) is the existence of high- momentum components for momenta 
q > 2 fm _1 due to the short-range and tensor NN correlations. This has been shown also for particular 
nuclei such as 2 H, 3 ' 4 He, 12 C, 16 0, 40 Ca and others and in nuclear matter within various theoretical 
correlation methods (see e.g. [f20|-[40"f and reviews in jt4|,[l5)). In addition, both experimental and the- 
oretical analyses confirm the conclusion that the high-momentum behaviour of the nucleon momentum 
distribution (n(q)/A at q > 2 fm _1 ) is similar for nuclei with mass numbers A— 2, 3, 4, 12 and 56 and 
for nuclear matter. 

The first aim of the present work is to apply the LFD method to calculations of the momentum 
distribution in the deuteron. Thus this study is an important complementary test of the possibility of 
the LFD to describe simultaneously both deuteron charge form factors (that has been shown in R]) and 



momentum distribution. In this line the results of calculations using wave functions of existing NN 
potentials are also given and compared with the LFD results and the YSD. 

The second aim of the work is to develop a method which would allow us to calculate the nucleon 
momentum distribution for a wide range of (A, .Z)-nuclei on the basis of the knowledge of n(q) in the 
deuteron and especially of its high-momentum components. The latter are important because in various 
processes with high transferred momenta only relativistic nucleons can take part ||. Such analyses 
require the knowledge of the wave function (WF) at momenta with magnitude of the nucleon mass. We 
make an attempt to relate the relativistic effects in n(q) in the deuteron with the explanation of the 
high-momentum components of n{q) in heavier nuclei which are important for calculations of the cross 
sections of various nuclear reactions. This question is discussed together with the effects of short-range 
NN correlations on the nucleon momentum distributions. 

A brief review of the LFD method is given in Sect. II. The results of the calculations of n{q) in deuteron 
within the LFD and the comparison with calculations using wave functions corresponding to existing NN 
potentials and the YSD are given in the same Section. In Sect. Ill a method to calculate the nucleon 
momentum distribution in (A, Z)-nuclei is developed using the natural orbital representation and the 
results on the high-momentum components of the nucleon momentum distribution in the deuteron. Sect. 
IV contains the conclusions of the present work. 



II. NUCLEON MOMENTUM DISTRIBUTION IN THE DEUTERON WITHIN THE LFD 

METHOD 



The results of the experimental and theoretical studies of the structure and form factors of the deuteron 



have been presented in the most recent review 41 . In it special attention has been paid to the different 
theoretical models for investigating the deuteron electromagnetic form factors and nucleon momentum 
distribution. Among them, the LFD method has been discussed as one of the promising approaches. We 
start this Section with a brief description of the LFD (e.g. ]7|- |l0| , |4^ -^4| ) which is applied to the analyses 
of relativistic bound systems at high relative momenta. The LFD is a self-consistent method in which the 
relativistic wave functions are the Fock column components of the state vector defined on an arbitrary 
light-front (LF) hypersurface given by invariant equation ujx = 0. The four-vector uj = (uq^lj), uj 2 = 0, 
determines the position of the LF surface. The state vector satisfies a dynamical equation following from 
rotations of the light front. In LFD the transformations of the reference system and of the light front 
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are independent of each other. The dynamical dependence of the state vector on the LF position is 
parametrized by the explicit dependence of the wave functions on the four- vector u>. The dependence on 
u, however, is not a property of the observable amplitude. The form factors depend on four-momentum 
transfer squared only and do not depend on u>, though they are related to the w-dependent wave functions 
JlO| . The relativistic deuteron WF on the LF ^(q,n) depends on two vector variables: i) the relative 
momentum q and ii) on a unit vector n along uj. Due to this, the WF is determined by six invariant 
functions instead of two [S- and Z?-waves) in the non-relativistic case. Each one of these functions 
depends on two scalar variables q and z — cos(q,n). In LFD these six functions are calculated within 
the relativistic one-boson-exchange model. The kernel in the dynamical equation for ^(q, ft) has been 
calculated flC[| using Lagrangians of interaction of nucleon with pseudoscalar (ir, rf)-, scalar (a, 6)- and 
vector (p, ui) mesons using the parameters of the Bonn model (4^]. Thus in the LFD calculations the 
deuteron structure is determined by the relativistic nucleon-meson dynamics, supposing that the nucleons 
in the deuteron interact by exchanging relativistic mesons [Q, without using potential approximations. 

The relations of the covariant LFD method to other relativistic approaches, such as the standard LFD 
fl46| , the Bethe-Salpeter formalism |47| and its three-dimensional reductions, the quasi-potential equations 
fl48| and the Gross equation method (4j| are given in j|. For instance, it is shown in |5(]J^] that after the 
projection of the Bethe-Salpeter amplitude on the light front, the six components of the LFD deuteron 
WF are expressed through integrals over the eight components of the deuteron Bethe-Salpeter amplitude. 
Provided the NN interaction is the same, these approaches incorporate by different methods the same 
relativistic dynamics. The WF's in LFD are the direct relativistic generalization of the non-relativistic 
ones in the sense that they are still the probability amplitudes. Therefore they can be used in the 
relativistic nuclear physics (e.g. JhJ). Another advantage of LFD is that due to the explicit relativistic 
covariance of the method, the form factors do not depend on the system of reference and on the direction 
of the z-axis. One should emphasize that the LFD WF's have been successfully used also in the hadron 
physics H. 



As mentioned in the Introduction, the LFD results give a good description of the data 41 for the 
tensor polarization £20 and for the structure function A(Q 2 ) (up to Q 2 — 3 (GeV/c) 2 ). Concerning the 
structure function B(Q 2 ) its minimum occurs below the position indicated by the data. The latter is 
related [0,^J to the insufficient accuracy of the perturbative approach in calculating the components of 
the LFD wave functions. We would like to emphasize, however, that the LFD method describes very well 
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the charge and quadrupole form factors which are expressed by the structure functions A(Q 2 ), B(Q 2 ) 
and the tensor polarization i 2 o- 

In this Section we firstly calculate the nucleon momentum distribution in deuteron n(q) within the LFD 
using the WF (^(q, ft)) which is normalized generalizing the non-relativistic normalization condition ]To| : 

1 m f _ T 2 _ d 3 q to f _,,_g _ -^d 3 q 



where F is expressed by the six scalar functions fi-e depending on the scalar variables q 1 and ft ■ q: 

F = fl + /| + fl + (3 + z 2 )/ 2 + (1 - ^ 2 )(/| + / 2 ) + (3z 2 - l)hh + 4z/ 4 (/ 2 + / 3 ) • (2) 



In (|l|) e(g) = \/ q 2 + to 2 , where to is the nucleon mass. As shown in |L0| , in the non-relativistic limit 
(q << to) the functions /3_e become negligible, /i^ do not depend on z and turn into S- and D-waves 
(fi~us, /2~~ u d) an< i the WF ^(q, ft) becomes the usual non-relativistic wave function. One of the most 
important properties of the functions /i-6 found in [ fiof is that for q > 2—2.5 fm -1 the component fs 
(being related mainly to 7r-exchange) exceeds sufficiently the S- and D- waves. This fact is very important 
in the calculations of n(q) in deuteron as it will be shown below. 

In JhJ a Legendre expansion of fis(q, z) relative to z has been used: 

/i(?,*)=E(2J + l)/i(?)flW (3) 
i 

(1=0,2 for i=l,2,3,5 and Z=l,3 for i=4,6) with given values of the coefficients f\{q). 

We calculate in our work the angle-averaged nucleon momentum distribution in the deuteron defined 
as: 



n(q) = C(q)F(q), (4) 

where 

C ^ = T^WT ( 5 ) 



and F(q) is the angle-averaged function F(q, z = cos( 



1 



F(q) = — J F(q,cos6)dn. (6) 
In accordance with Eq. (|l|) the normalization of n(q) is given by: 

n(q)d 3 q = l. (7) 



The LFD calculations have shown that, as expected, the most important contributions to the total 
n(q) give terms related to the fx, f 2 and f 5 functions 

n(q)~nx(q) +n 2 (q) +n 5 (q), (8) 

where 

nx(q)=C(q)JW), n 2 (q) = C(q)J$(qj, n 5 (q) = C{q){\ - z 2 )/|(g). (9) 

Here we would like to note that the use of the functions fx, f% and f§ averaged over z = cos{q, ft) can be 
justified by their smooth (almost constant) behaviour as functions of z at different values of q (see Figs. 
10 and 11 of |fjofl ). The contributions of nx, n 2 , nx2 = nx + n>2 and ng are compared in Fig. 1. It can be 
seen that, while the functions fx and f 2 give a good description of the YSD of n(q) for q < 2 fm _1 (like 
the S- and £>-wave functions in the non-relativistic case) , it is impossible to explain the high- momentum 
components of n(q) at q > 2 fm -1 without the contribution of the function fy. We note that the deviation 
of the total n(q) from the sum nx2 = n\ + ri2 starts at q around 1.8 fm _1 . All this shows the important 
role of NN interactions which incorporate exchange of relativistic mesons in the case of the deuteron. We 
emphasize the consistency of the LFD method obtained in the simultaneous description of both deuteron 
momentum distribution (made in this work) and of charge form factors jjj. Here we would like to discuss 
this point. Considering the LFD results on the nucleon momentum distribution in the deuteron we do 
not imply that the LFD method and the contribution of the f$ function are the only way to understand 
the high-momentum components of n(q). In the following of this Section we will pay attention to the 
results of other approaches considering, however, this question in close connection to the possibility of 
simultaneous description of both momentum distribution and form factors. It is known, for instance, 
that the deuteron wave function related to the Argonne Vxs potential for the NN interaction |l2] gives a 
realistic description of n(q) in the deuteron ]5^| . However, as can be seen from Figs. | and | of Ref. § , 
the recent non-relativistic impulse approximation calculations using the Argonne vig potential and even 
those including meson exchange currents and relativistic corrections do not give a good agreement with 
the new data of Q on the tensor polarization t 2 o and the charge form factors Gc and Gq. 

We calculate in the present work also the angle-averaged nucleon momentum distribution in the 
deuteron using S- and D-wave functions ^s{q) and ^d(q) corresponding to various NN potentials, 
such as the charge-dependent Bonn potential | p3[ , the Argonne Vxs [0, the Nijmcgen - I ,- II and - Reid 
93 j54j and Paris 1980 |55| in the expression: 
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n (i) = ^: [*!(«) + = + nu(g) (10) 

with 

Jn(q)d 3 q=l. (11) 

In Fig. U the result for 71(5) using the charge-dependent Bonn potential |53j | is given and compared with 
the YSD. As can be seen, the D-component of n(q) is important but even its inclusion does not give a 
very good agreement with the data for q > 2 fm _1 . In the next Fig. || we present the LFD result for 
n(q) compared with the calculations using the WF's corresponding to Nijmegen-I,-II, -Reid 93, Argonne 
vis and Paris 1980 NN potentials and with the y-scaling data. We would like to note that: i) the results 
of the calculations using the NN potentials, such as Nijmegen-II, -Reid 93, Argonne v%s and Paris 1980 
(shown in Fig. [5j) are in better agreement with the YSD than those using the charge-dependent Bonn 
potential (Fig. ||) . This might be related to the fact that these potentials describe NN phase shifts up to 
larger energies (e.g. the Nijmegen-II potential gives reasonable pp phase shifts up to 1.2 GeV, while the 
charge-dependent one-boson exchange Bonn potential fits the phase-shift data below 350 MeV); ii) It can 
be seen from Fig. || that there are small differences between the curves corresponding to different NN 
potentials for q < 3 fm -1 (which give a good description of the YSD and almost coincide with the LFD 
result) and larger ones for q > 3 fm -1 . Large differences take place, however, between all of them and 
the LFD result for q > 3 fm . An important conclusion of this fact is that data at higher momentum 
transfers are needed (and probably a more refined analysis of the currents used for the coupling of the 
electron with the nucleons) in order to distinguish the properties of the covariant LFD method from 
the potential approaches. The momentum region of q just larger than 3 fm -1 will be decisive in this 
sense. The main advantage of the LFD method is, however, as mentioned already, that it describes 
simultaneously the deuteron charge form factors as well, which is not the case with e.g. the Argonne vis 
NN potential method (see Fig. 4 in jy) . 

III. NUCLEON MOMENTUM DISTRIBUTIONS IN 4 He, 12 C AND 56 Fe 



As mentioned in the Introduction, realistic many-body calculations (e.g. |20|- f40| , |l4|l5{ ) have shown the 
existence of high-momentum components in the nucleon momentum distributions for studied nuclei at 
q > 2 fm -1 , due to presence of short-range and tensor correlations. This conclusion agrees well with 
the y-scaling data |l6]-[l^] . Secondly, it has been shown that all nuclear momentum distributions ha (<?) 



for (A, Z)-nuclei (A > 2) at these high momenta are simply rescaled versions of the nucleon momentum 
distribution in the deuteron n(q) ]56|] : 

n A (q)=a A n{q), (12) 

where a a is a constant. 

Since the magnitude of the high- momentum tail of n A (q) is proportional to the mass number, this 
effect is associated with the nuclear interior rather than with the nuclear surface. In the present paper 
we develop a model to calculate n A (q) for nuclei with A > 2 from that of the deuteron. The method has 
similarities to that suggested in [^7j on the basis of 4 He. However, in the present work we incorporate 
the analysis performed in Sect. II for the deuteron as the simplest two- nucleon bound system to the 
consideration of riA(q) for heavier nuclei, namely because the high- momentum components are related 
to the specific features of the two-nucleon interaction at short distance (around the repulsive core). The 
latter determine the corresponding behaviour of the nucleons in the central part of the nucleus, where the 
density is higher. There the two nucleons could be closer to one another and the short-range correlations 
could be operative. In this sense, we consider the present method as a next step and improvement on 
that of |7). 

In this Section we suggest a practical method to calculate the nucleon momentum distribution for the 
general case of (A, Z)-nuclei (^4 > 2), while as we mentioned already, the correlation methods consider 
only particular nuclei (e.g. 3,4 He, 12 C, 16 0, 40 Ca and nuclear matter) and due to their complexity the 
applications to other nuclei is often a very difficult task. 

We use the natural orbital representation (NOR) |5|,[l4|,[l5| of the one-body density matrix. In our 
case we consider its diagonal elements in the momentum space, which represent the nucleon momentum 
distribution n A (q). The NOR allows us to use the transparency of the single-particle picture even being 
beyond the mean-field approximation and to account for the correlation effects in a proper way. Using the 
common theoretical ground of NOR, the method makes it possible to combine the mean-field predictions 
for riA(q) (which are expected to be realistic at q < 2 fm _1 ) with that part of the momentum distribution 
which includes correlation and relativistic effects. 

In the natural orbital representation the momentum distribution in a nucleus with A nucleons can be 
written in the form: 

n A (q)=N A [n h A (q) + n p A (q)}, (13) 
where N A is the normalization constant and 
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F.L. 

n h A (q) = jJ2 2 ( 2 J + l)KijC(q)\R nlj (q)\ 2 (14) 

nlj 



oo 

= ^ E 2 ( 2 -?' + l)A„yC( ? )|fi n y(«)| 2 (15) 

F.i. 



are the hole- and particle-state contributions, respectively. In ( pL4[ ) and (15) F.L. denotes the Fermi level, 
C(q) is given by (|^), A„;j are the natural occupation numbers (0 < X n ij < 1) for a state with quantum 
numbers nlj and R n ij{q) are the natural orbitals. We call hole-state natural orbitals those natural orbitals 
for which the numbers X n ij are significantly larger than the remaining ones, called particle-state natural 
orbitals ]5S[ |. As shown in p9| , the high- momentum components of the total n A (q) caused by short-range 
correlations are almost completely determined by the contributions of the particle-state natural orbitals. 
This fact, together with the approximate equality of the high- momentum tails of n(q)/A for all nuclei 
allows us to make the main assumption of the method that the particle-state contributions to n A (<?) are 
almost equal for all nuclei. The deuteron n(q) from (^) can be written in a form similar to (|l3|): 

n(q)=N d [n h d (q) + n p d (q)} (16) 

with 

n h M = \[ni{q)+n 2 {q)l (17) 



n P M = ^n 5 (q), (18) 



and N d = 2. The assumed equality of the particle-state contributions for all nuclei implies the equality 



of the particle-state contribution n p A (q) for (^4, Z)-nuclci from (15) with that to the deuteron momentum 
distribution n p d (q) from ( |l8| ) 

n p A {q)^n%q). (19) 

The right-hand side of ( [l9] ) can be taken from different sources. This could be the YSD or results 
of calculations within given theoretical correlation models. Having the results from Sect. II we take 



n p d (q) to be related to n$(q) (18) because within the LFD namely this function is responsible for the 



high-momentum tail of n(q) for deuteron: 

n p A (q)^n 5 (q). (20) 
Then the resulting momentum distribution for (A, Z)-nucleus normalized to unity has the form: 



F.L. . 
n A (q) = N A [J2 2(2j + l)Xnl 3 C(q)\R nl] (q)\ 2 + jMq)} , (21) 

nlj 

where 

F.L. 

N A = {4tt / dq-q 2 ^ 2(2j + l)A ni ,C(g)|i?„ (g)| 2 + -^(g)]}" 1 ■ (22) 

It is known from that the hole-state natural orbitals are almost unaffected by the short-range cor- 
relations and, therefore, the functions R n ij(q) in ( piT| ) and (|2^) can be replaced by single-particle wave 
functions from the mean-field approximation (shell- model or Hartree-Fock WF). In our calculations we 
use Woods-Saxon s.p. wave functions for protons and neutrons. The hole-state occupation numbers X n ij 
are close to unity and we set them equal to unity with good approximation. These properties lead to a 
similarity of our model to that suggested for calculations of the spectral function . 

The calculations of n A (q) for 4 He, 12 C and 56 Fe are given in Figs. Q-^ and are compared with the 
y-scaling data from |lq| . We note that, as expected from (|l|) the value of a A must be almost constant. 
In our case a A — N A ■ A/2 and its numerical values are 8.335, 8.352 and 8.372 for 4 He, 12 C and 56 Fe, 
respectively. 

The dashed line in Figs. |||| gives the calculated normalized hole-state contribution (the mean-field 
component) to n A (q) only, i.e. without the contribution of n^{q) in Eq. (pl|). The comparison between the 
total n A (q), its hole-state contribution and the y-scaling data clearly shows the role of ns(g) contribution. 
It can be seen that the difference between the results with (q) and without it starts for all three nuclei 
at q ~ 1.5 fm" 1 and increase rapidly at q > 1.5 fm . This shows the role of the high- momentum 
components of the deuteron momentum distribution 715 which can be considered as an "ingredient" of 
the high- momentum components of the momentum distribution of (A, Z)-nucleus. In our opinion, the 
related function /s of the LFD method incorporates the main part of the short-range features of the 
NN interactions which determine the correlation effects seen in n A (q) when calculated within the non- 
relativistic correlation methods. 

The good agreement of the results on n A {q) with the data for the three nuclei achieved, including 
also the region of the high momenta, proves the assumption that the nucleon momentum distribution 
in nuclei can be extracted on the basis of the realistically described properties of the deuteron as the 
smallest bound system of two nucleons. 
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IV. CONCLUSIONS 



The results of the present work can be summarized as follows: 

(i) The light-front dynamics method is applied to calculate the nucleon momentum distribution in the 
deuteron. A very good description of the y-scaling data is obtained. It is shown that, while the functions 
/i and fi describe well the YSD of n(q) for q < 2 fm _1 , it is impossible to explain the high- momentum 
components of n(q) at q > 2 fm _1 without the contribution of the function f$. Thus, this study gives 
a positive answer to the question about the possibilities of the LFD method to describe simultaneously 
both deuteron momentum distribution and charge form factors (the good agreement with the latter as 
well as with the data for the tensor polarization £20 and for the structure function A(Q 2 ) up to Q 2 ~ 3 
(GeV/c) 2 was shown in @Q). 

(ii) The nucleon momentum distribution in deuteron was calculated also in this work using the S- and 
D-wave functions corresponding to various NN potentials (charge-dependent Bonn, Paris 1980, Nijmegen- 
I and -II, -Reid 93, Argonne vig). We show that the results using these potentials explain very well (with 
the exception of the charge-dependent Bonn potential) all the available data for n(q) up to q ~ 3 fm -1 
exactly like the LFD method. However, as shown in jl| , the recent non-relativistic impulse approximation 
calculations (e.g. with Argonne vis potential) and even those including meson exchange currents and 
relativistic corrections do not describe well the new data Q| on the tensor polarization £20 and the charge 
form factors Gc and Gq. This is the main difference between the results of the LFD method and of 
the calculations using WF's corresponding to NN potentials for the important deutron characteristics 
considered. Secondly, we should mention also that for q > 3 fm -1 the LFD results for n(q) deviate 
strongly from those of the calculations using NN potentials. Thus data at higher q are needed in order to 
distinguish the abilities of the LFD from those of the potential methods. Here we would like to mention 
also other relativistic methods such as the method in [p"l| (using a three-dimensional reduction of the 
Bethe- Salpeter equation) and the method in |n| (based on the Gross equation |l9|]) which could give a 
simultaneous description of both deuteron form factors and momentum distribution. 

(iii) A correlation method for calculating the nucleon momentum distribution in nuclei with A > 2 
nA(q) is proposed. The method uses the natural orbital representation of the one-body density matrix 
and combines the mean-field component of riA{q) with the correlated high-momentum components taken 
from the deuteron momentum distribution. The method is based on the well-known results of realistic 
many-body calculations which have shown that, due to the short-range and tensor NN correlations, 
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all nuclear momentum distributions riA^q) at high momenta {q > 2 fm -1 ) are rescaled versions of the 
nucleon momentum distribution in the deuteron n(q). In the present work we incorporate in the method 
the LFD results for n(q). The role of the wave function f$ is clearly shown by the comparison of the 
calculations of nA{q) using the hole-state (mean- field) contribution without and with an inclusion of 
the 715(75 ) contribution. In our opinion, the LFD function f$ incorporates the main part of the short- 
range features of the NN interactions which can be seen in calculations of nA^q) within non-relativistic 
correlation methods. The suggested method gives an easy practical way for calculations of nA{q) for any 
nucleus because only hole-state wave functions and occupation probabilities enter its main relationships 
(Eqs.(|2l|) and (p2|)). The method gives a good description of the y-scaling data for nA{q) in 4 He, 12 C and 
56 Fe including the high-momentum region, showing a clear difference from the mean-field predictions for 
q > 1.5 fm -1 . 
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FIG. 1. The nucleon momentum distribution in deuteron. The contributions of ni, n^, = ni + and n$ 
are presented. The {/-scaling data are from The normalization is: J n(q)d?q = 1. 

FIG. 2. The nucleon momentum distribution in deuteron (solid line) calculated using Eqs. © and (0) with 
S- and D-wave functions corresponding to the charge-dependent Bonn potential |k| . S- and D-contributions are 
given by dashed and dotted line, respectively. The {/-scaling data are from |l8| . 

FIG. 3. The nucleon momentum distribution in deuteron calculated with the LFD (solid line) in comparison 
with the results presented by numerated arrows, as follows: 1-Argonne vis, 2-Nijmegen Reid 93, 3-Nijmegen I, 
4-Nijmegen II, 5-Paris 1980 NN potentials and with the {/-scaling data [ p^[ . 

FIG. 4. The nucleon momentum distribution in 4 He calculated using Eqs. (pll) and (||) (solid line). The 
dotted line represents the hole-state contribution only. The {/-scaling data are from ifjjjfl . The normalization is: 
/ n 4 (g)d 3 g = 1. 



FIG. 5. The same as in Fig^ for 12 C. 



FIG. 6. The same as in Fig.fi] for 56 Fe. 
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